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Abstract 

Mass parameters for collective variables of a dinuclear system and strongly de- 
formed mononucleus are microscopically formulated with the linear response 
theory making use of the width of single particle states and the fluctuation- 
dissipation theorem. For the relative motion of the nuclei and for the degree 
of freedom describing the neck between the nuclei, we calculate mass param- 
eters with basis states of the adiabatic and diabatic two-center shell model. 
Microscopical mass parameters are found larger than the ones obtained with 
the hydrodynamical model and give a strong hindrance for a melting of the 
dinuclear system along the internuclear distance into a compound system. 
Therefore, the dinuclear system lives a long time enough comparable to the 
reaction time for fusion by nucleon transfer. Consequences of this effect for 
the complete fusion process are discussed. 

PACS:25.70.Jj, 24.10.-i, 24.60.-k 
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I. INTRODUCTION 



The synthesis of transuranium and superheavy elements, the production of superde- 
formed nuclei and nuclei far from the line of stability stimulate the study of phenomena of 
complete fusion in heavy ion collisions at energies less than 15 MeV/nucleon. The fusion 
of heavy ions is characterized by the formation of a dinuclear system (DNS) at the initial 
stage of reaction when a significant part of the kinetic energy is transferred into excitation 
energy. In that case one has to describe the evolution of the DNS and strongly deformed 
nucleus by collective coordinates Q, namely by the distance R between the centers of collid- 
ing nuclei (or relative elongation A of the system), the mass and charge asymmetry degrees 
of freedom rj = {Ai — A2)/A and rjz = {Zi — Z2)/Z, respectively, {Ai, A2 and Zi, Z2 are 
mass and charge numbers of the nuclei, A = Ai + A2, Z = Zi + Z2), the parameter e of the 
neck and deformation coordinates of nuclei Then dissipative, conservative and inertial 

forces for collective variables, arising from the nucleus-nucleus interaction, are needed to be 
determined. 

The present fusion models can be distinguished by the choice of the relevant collective 
variables along which fusion mainly occurs. While many models study the fusion in R at 
a practically fixed value of rj, the DNS model |0,|lO| considers the DNS evolution in mass 
asymmetry by nucleon or cluster transfers as the main path to the compound nucleus. The 
DNS model assumes basically that the neck degree of freedom is fixed in the evolution in rj 
and the nuclei are hindered to melt together by a variation in the relative distance. As was 



shown in Refs. |[TT|-[14| this occurs due to the relatively large inertia of the neck degree of 



freedom and structural forbiddenness effects. In this paper we study dynamical restrictions 
for the growth of the neck in the DNS and suggest proper methods of the calculation of the 
DNS inertia tensor. 

There are various macroscopical and microscopical approaches to calculate the inertia 
tensor |T^. The macroscopical approaches (see, for example, PJ^) are based on the hy- 
drodynamical model of the nucleus. A calculation of the inertia tensor with a theory for 



quantum fluid dynamics is suggested in Ref. [16|. By using a random- matrix model to de 



scribe the coupling between a collective nuclear variable and intrinsic degrees of freedom and 



with the help of the functional integral approach, mass parameters are derived in Ref. [jTT 
In the linear response theory []TB|,|T^ the inertia tensor is found for fissioning nuclei. The 
microscopical approaches mainly use the cranking type expression and perform calculations 
in different single particle bases applying adiabatic |l3|J^,pO|] or diabatic |l2T| , ^ two-center 
shell models. Difficulties in the cranking type calculations arise for collective motions with 
large amplitudes, for example, in fusion or fission, due to pseudo-crossings or crossings of 
levels in the single particle spectrum. Some publications disregard the contributions from 
the crossings (pseudo-crossings) which means a neglect of effects of configuration changes on 
the mass parameters during the evolution of the nuclear shape in spite of the fact that the 



collective inertia is strongly influenced by level crossings (pseudo-crossings) [0,0] . In order 



to overcome this problem, two-body collisions should be regarded which lead to a width of 
the single particle levels and an effective reduction of the level crossing effects. For example, 
calculations of the nuclear inertia in a generalized cranking model with pairing correlations 
yielded masses of about one order of magnitude larger than the ones without pairing [p5| . 
One of the aims of the present paper is to obtain analytical expressions for mass parame- 



2 



ters using models and methods which include residual interaction effects. In Sect. 2 the mass 
parameters are obtained within the linear response theory taking the ffuctuation-dissipation 
theorem and the width of single particle states into account. The same mass parameters 
are also derived by Fermi's golden rule and by smoothing the single particle spectrum in 
the mean-field cranking formula. In Sect. 3 the mass parameters for the relevant collective 
variables (mass asymmetry, elongation, neck and deformation parameters) of the DNS and 
strongly deformed nuclear systems are evaluated in the two-center shell model with adiabatic 
and diabatic bases. Sect. 4 contains a short summary. 



II. MICROSCOPICAL INERTIA 

A. Derivation from collective response function 

Let us consider a nuclear system described by a single collective coordinate Q and intrinsic 
single particle coordinates Xi (with the conjugated momentum pi) and assume the following 
effective Hamiltonian 



H{xi,pi, Q) = H{xi,pi, Qo) + {Q - Qo)F{xi, pi, Qo) 

1.^ .2 ^ d^H{x,,pi,Q) 
+ 2 w - Woj < >Qo,ro (ij 

The shape of the nuclear mean field is changed with the collective coordinate Q that intro- 
duces the coupling between Q and the nucleonic degrees of freedom. Eq. (1) is obtained by 
expanding the Hamiltonian to second order in the vicinity of Qo and by approximating the 
second order term in a kind of unperturbed limit. The coupling term between the collective 
and intrinsic motion is proportional to the first order in 6Q = Q — Qo with an operator F 
given by the derivative of the mean field with respect to Q in the neighborhood of Qq. In 
this way global motion is described within a locally harmonic approximation. 

By applying the linear response theory |1T8| , ^ , the Fourier transform of the collective 
response function is defined as 

where x{^) is the Fourier transform of the response function for intrinsic motion and mea- 
sures how, at some given Qq and temperature Tq, the nucleonic degrees of freedom react to 
the coupling term FSQ. The coupling constant k is written in the form 

, d^H{x„Pi,Q) d^E{Q,So). , , 
-k =< >Qo,To= Iqo + x(^ = 0) = C(0) + x(0) (3) 

with x(0) and C(0) being the static response and stiffness, respectively. Since the constant k 
is entirely determined by quasi-static properties, it is no surprise that E is the internal energy 
at a given entropy 5*0 or the free energy at a given temperature Tq. The structure of Eq.(2) 
reflects self-consistency between the treatment of collective and microscopic dynamics. It 
expresses the response of the system of interacting nucleons in terms of the response of the 
individual nucleons. 
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The local motion in the Q variable is described in terms of the collective response function 
Xcoii- A short calculation leads to the following approximation for the mass coefficient 
g|T9|^,|23: 



M 



1 d\xcou{uj))-\ 



2P 



1 1^=0 



C^(0),2r,.c. , 7^(0)- 



x(o) 



X(0) 



where 



1 d\{u) _ 



(4) 



(5) 



is the inertia in the zero-frequency limit of the second derivative of the intrinsic response 
function. M'^'' can be shown to be similar to the one of the cranking model. For many 
applications, the value of C(0)/x(0) is much less than unity. The additional term 7^(0)/x(0) 
in Eq. (4) gives a positive contribution to M where 7(0) is the friction coefficient defined by 



7(0) 



.dxiuj) 



dx\uj) 



2Tr 



■/(o). 



(6) 



The dissipative part of the response function x" i^) is connected with the dissipative part 
of the correlation function ijj (a;) through the fluctuation-dissipation theorem: 



X [u}) = -tanh{—)tfj [uj). 



(7) 



The il>" (uj) has a singularity of (5-function type at u = 0: 

^"M = 2vr^°(5M+V^^H (8) 
with being regular at u = 0. In the case of an independent particle model we have 



■n{ej)[l 



n{ek)][6{huj - ekj) + S{huj + e^j)]. 



(9) 



Here, Ckj = Ck — Cj is the difference of single particle energies, n{ej) are the occupation 
numbers and Fjk =< j\F\k > the single particle matrix elements of the operator F. At 
j = k and u = we flnd the contributions from the diagonal matrix elements: 



= E \Fkk\Mek)[l - n{ek)] = Tq^ 

k k 



dn{e) 



de 




(10) 



The last part in (|TOD was derived with a Fermi distribution for the occupation numbers, 
which is characterized by the temperature Tq. The value of Tq does not effectively go to zero 
with decreasing excitation energy because each single particle level has a width due to the 
two-body interaction. Indeed at zero excitation energy the distribution of the occupation 
numbers deviates from a step function at least due to pairing correlations. If we replace the 
5-functions in Eq.(8) or (9), we have to apply a Lorentzian functions with the double single 
particle width 2r where F is the single particle width, because hu is the transition energy 
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between two single particle states [|T8|. Therefore, we substitute the 5— functions in Eq. (9) 
by T /['K{{huo ± CkjY + r^)]. Then using Eqs. (6)-(10), we can write the friction coefficient in 
the following form: 



nondiag / 



where 



h 



dn{e) 



de 



e=efe 



dQ. 



(11) 



(12) 



For smaller temperatures Tq < 2 MeV which are of interest here, ^'^^'^^(^q^ jg much larger 
than 7'^°"-'^*'^5(^o) |]TB[. The static response is found as 



X(0) = hm / = hm / ^ ^..a.(o) + x-'-m, (13) 



J Ti ^ — e^o J hn uo — ie 

-oo — oo 



where 



/-3(0) = ^ 



(9n(e) 



de 



e=efe 




(14) 



With realistic assumptions 7'='^'^^'(0) > 7"°'^'^^«9(0) and x'^'"^(0) > ;;^"°"'^^'^9(0) and neglecting 
C(0)/x(0), we can divide the mass parameter (4) as 



_ jy^diag _|_ j^nondiag 

The contribution of the diagonal matrix elements of F to M are 

( dekV 



dQ 



(15) 



(16) 



If the single particle widths are properly taken into account, the nondiagonal contributions 
to the inertia are 12011 



j^nondiag _ j^cr _ ^2 



k^k' '^kk 



, + r2 



(17) 



The main contribution to M is the diagonal part Af^*"^ because it dominates for collective 



variables which are responsible for changes of the nuclear shape of the system p3|,^,^ . Note 
that the calculation of M'^'^"'^ is simpler than M^°^'^^"-^ . For the case that the pairing residual 
interaction is regarded and only diagonal matrix elements in the cranking formula are taken 
into account, Eq. (p!6| ) was obtained with F = A (A is the pairing gap) in Ref. p|,p5|. Start- 
ing with an equation for the single particle density matrix extended with an approximate 



incorporation of particle collisions in the relaxation time approach, the authors of Ref. 
derived an expression similar to ([T6|) (with F = /i/r, r is the relaxation time) but with a 
negative sign. This negative sign arises from the fact that the condition of self-consistency 
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between collective and nucleonic dynamics was disregarded in which is important for a 
correct calculation of the mass parameters. It was stressed in [l^,!^] that within the linear 
response theory the diagonal component of the friction parameter originates from the "heat 
pole" of the correlation function and vanishes when the system is ergodic. As shown in 
Ref . |31[] , the well necked DNS-type configurations are not ergodic and stable against chaos. 
Even at zero excitation energy the level crossings at the Fermi surface lead to consider- 
able mass flow |]23|J2^,|2S| and the diagonal component of the correlation function (or mass 
parameter) does not vanish. 

Besides the mass and friction coefficients, the diffusion coefficients D^i {k,l = {Q,P)) 
must also have a component diagonal in the matrix elements of F because they are connected 
with correlation functions. For example, the diffusion coefficient in momentum is defined 
as: 



D 



pp 



1 i"( 

-w (UJ 



o) = ro7(o). 



B. Derivation from Fermi's golden rule 



By setting Hj = {Q — QQ)F{xi,Pi, Qq) as perturbation (see Eq. (1)), the decay rate of a 
collective state \n > with energy En to the collective state \m > with energy is given in 
lowest order according to Fermi's golden rule: 



27r 

w(n m + hu) = —I < m\Q — Qo\n > P 
n 



X / d{huj)\ < fiuj\F\0 > \ 6{E„, - Em - huj)pgs{fiuj). 



(18) 



Here, the integral is taken over the final states of the intrinsic system with the density Pgg. 
With the fluctuation-dissipation theorem for small temperatures we have p2| : 



< huj\F\0 > \^ pqs{hu)d{?iuj) = — —R{uj)duj 



(19) 



with the relaxation function R{uj) = x"{^)/^- The half-decay width is obtained from 
Eq.(|l8D as 



(20) 



Using the properties of the response function and a Taylor expansion 



R{uj) 



x"(^ = o) + 



duj 



2 duj^ 



uj=0 



(21) 



and the standard formula for mass M {En 7^ Em) 



M = ^(E \<m\Q- Qo\n > \^[En - Em])'' 



(22) 
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we obtain by setting F = r„ with Eqs. (IT 



h dy" (uj) , h , . , , 



Large temperatures in (19) effectively lead to a temperature dependence of T in (P^D- Since 
7(0) in Eq. (11) contains the terms with the diagonal matrix elements of the operator F, the 
mass parameter M also has the diagonal component M'^*"^ (HID- So, the contributions to 
the mass parameter can be again classified as those with diagonal and nondiagonal matrix 
elements, respectively. 

That the mass parameter is proportional to the friction coefficient (see Eq. (0)), has an 
analogy in the hydrodynamic model. For multipole moments A of the nucleus with A > 1, 



the following ratio in the limit of an irrotational flow was derived in [SB 

3^2/3 1 



irr 



^ '^^ 87r(2A + l)(A-l)roz/' 
where p is the coefficient of the two-body viscosity and tq = 1.2 fm. 



C. Derivation from the mean field cranking formula 

Using the single particle spectrum and the corresponding wave functions, one can obtain 
the mass parameter with the cranking formula 

In reality the Hamiltonian of the system contains a residual two-body interaction between 
the nucleons in addition to the mean field. The residual coupling distributes the strength 
of single particle states over more complicated states. This spectral smoothing causes the 
effect that the sum over a and (3 appearing in (|2^) also includes diagonal terms with a = (3. 
Let us prove this statement. 

The Eq. (|2^) can be rewritten as 

d |2"-(ei) - n{e2) 

oQ 62 - 1\ 



= E / ^ei5(ei - e„) / de2K^2 - < e,\^\e2 > plH^lLJlM. (25) 

^ J J 0(J 69—61 

Next we use the following replacements 

J deig{ei) E' / de2gie2) E' (26) 



Sie - 6.) - p,(6) = ^ (^_^^).^(r/2)2 (27) 



and the approximation 
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D^\<h\-^\k2>\^ ^\Fu,k,?. (28) 

Here, D = 1/g is average energy distance between single particle states. Then we express 
the mass (B5|) as 



The energy spreading of the single particle states is taken into account in Eg. ([29|) . A line 
broadening happens if collisions of particles and holes with the background result in single 
particle and single hole strength functions that are concentrated around the original single 
particle energies. The quantity Dpk determines the average strength function for a particle 
in state k pi]. In the hmit e^i = = Cfc, for Ca^/s ~ e^, i.e. when two neighboring levels 



near the level k are considered, and S/vr^ ^ 1, Eq.(|^) leads to Eq. ([T6|) . Thus, diagonal 
terms in the mass parameters appear because of the finite width of the single particle levels 
due to the residual interaction. 



III. RESULTS OF CALCULATIONS 

A. Adiabatic two-center shell model 

Since in fusion and quasifission we deal with strongly elongated systems, the two-center 
shell model (TCSM) [^^,|rT| is appropriate for calculating the potential energy surface. In 
the TCSM the nuclear shapes are defined by the following set of the coordinates: The 
elongation A = 1/{2Rq) measuring the length / of the system in units of the diameter 2Rq 
of the spherical compound nucleus and used to describe the relative motion, the mass and 
charge asymmetry coordinates rj and rjz-, respectively, the neck parameter e = Eq/E' defined 
by the ratio of the actual barrier height Eq to the barrier height E' of the two-center oscillator, 
and the deformation parameters /3j = Oj/fej, z = 1, 2, of axially symmetric fragments, defined 
by the ratio of the semiaxes of the fragments. Since collisions above the Coulomb barrier are 
discussed, we firstly consider spherical nuclei with f3i = 1 and then analyze the deformation 
effects. 



In order to calculate the width of the single particle states, we use the expression |1T8|J35 

r„ 1 + |(es - )V (irT„)2]/c2 ' 
Here, ep is the Fermi energy. Both parameters Fq and c are known from experience with 



the optical model potential and the effective masses Their values are in the following 
ranges: 0.030 MeV"^ < Fq"^ < 0.061 Me 15MeV< c < 30MeV. For small excitations, 
Eq. (PPD is reduced to the expression known in the theory of a Fermi liquid. Since each 



single particle state has its own width, Eq. (|T^) is generalized as {Qi = X^rj^rjz, e, (3i, P2)- 

Mr^ft^ElH^. (31) 
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For the Fermi occupation numbers n{ek), the function 



duk 1 , _2 f ek - ep 



^ _ ^ ^ ^^^^ 

has a bell-hke shape with a width Tq and is peaked at the Fermi energy ep. 

Various calculations of the mass parameter for the motion in A were carried out with 



expressions similar to Eq. (JST]), for example in [1,20,25]. When the system adiabatically 
moves towards the compound nucleus, the value of Mw approximately increases by a factor 
10-15 in our and other calculations. In this paper we concentrate on the calculation of the 
mass parameter Mg^ for the motion of the neck to test whether the DNS exists sufficient time 
enough with a relatively small neck. The dependence of M^^ on e is presented in Fig. 1 for 
the system ^^°Pd+^^°Pd at A = 1.6 which corresponds to the touching configuration in this 
symmetric reaction. The obtained values of M^^ have the same order of magnitude as in [20] 
where the pairing correlations were taken into consideration. The value of M^^ increases by 
a factor 2.5 when the system falls into the fission-type valley [|lT|. This increase refiects the 



decrease of the shell correction SU with e towards e — > 0. Smaller values of 5U correspond 
to larger masses. 

In order to obtain the same TCSM potential energy as in the DNS model for the touching 



configuration, the neck parameter e should be set about 0.75 [11|. With this value of e the 
neck radius and the distance between the centers of the nuclei are approximately equal to 
the corresponding quantities in the DNS. 

For the parameter c in Eq. (^) we use the "standard" value 20 MeV since the masses 
(31) depend only weakly on this parameter. Setting the parameter Fg^ = 0.045 MeV~^ in 
(pT|) and comparing our results with obtained in the Werner- Wheeler approximation, 

we find = M^^^, M,, ^ (20 - 30)Mj^, M^, ^ 0.4Ml^^ and M^J^/M^jMTe < 1, 
practically independent of the mass number of the system. Therefore, we can conclude 
that the microscopical mass parameter of the neck is much larger than the one in the 
Werner- Wheeler approximation and the nondiagonal component M^^ is small. Since at the 
touching configuration the slope of the single-particle levels is small and changes slowly 
with decreasing elongation, the microscopical mass parameter in A is close to its smooth, 
hydrodynamical value. In contrast, a large amount of internal reorganization, which occurs 
at level crossings with decreasing e, leads to a large neck inertia of the initial DNS. So, 
the value of M^^ exceeds the mass in the hydrodynamical model due to large values of 
\dek/de\. The restriction for the growth of the neck may be understood by analysing the 



single particle spectrum as a function of e |TT[. For large e and well necked-in shapes, the 



single particle spectrum shows a good shell structure. The levels show an increasing number 
of level crossings with increasing e. 

The time-dependence of the neck parameter calculated with the microscopical and 
Werner- Wheeler mass formulas are compared in the upper part of Fig. 2. The lower part 
of Fig. 2 shows trajectories in the {e, A)-plane, calculated with microscopic and Werner- 
Wheeler masses. An adiabatic potential energy surface is used in all these calculations. 
Since there are no suitable barriers at smaller values of A and e in the adiabatic potential 
which hinder a growth of the neck, the neck parameter and system length decrease steadily 
to smaller values, faster in the case with the Werner- Wheeler masses and much slower with 
the microscopical masses. The experiments on fusion of heavy nuclei can not be explained 
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as a melting with increasing neck together with a decreasing A in an adiabatic potential 



rT| . It seems that there is an intermediate situation between the adiabatic and diabatic 
limits. The study of the transition between diabatic and adiabatic regimes gives a potential 
energy surface which contains quite high barriers for the motion to smaller A and e [|r3| , |36| . 
Therefore, the dynamical calculations with the adiabatic potential energy show a maximal 
possible growth of the neck. Due to the large moment of inertia in the heavy nuclear system, 
the dependence of potential energy on the angular momentum in dynamical calculations is 
disregarded here. Indeed, in fusion reactions with massive nuclei only low angular momenta 
(< 20 ^) contribute to the evaporation residue cross section. For the study of fusion, the 
isotopic composition of the nuclei forming the DNS is chosen with the condition of a N/Z- 
equilibrium in the system with rjz following the value of r] [10]. In processes developing in 
the shorter times, the coordinate rjz could be considered as independent collective variable. 

As result of Fig. 2 we find that the microscopical mass parameters keep the system near 
the entrance configuration for a sufficient long time comparable with the time of reaction 
even in an adiabatic potential. Therefore, this situation justifies the assumption of a fixed 
neck as we assume it in the DNS model [0,|l3]- When the DNS configuration exists a long 
time, then thermal fiuctuations in the mass asymmetry coordinate play the essential role in 
the fusion process. Indeed these fiuctuations are responsible for the fusion in the DNS model 
||7,10]. Thus, the dynamical restriction for the growth of the neck can be caused partly by 
a large microscopical mass parameter for the neck motion, partly by the potential energy 
surface intermediate between diabatic and adiabatic limits. 

In the adiabatic two-center shell model the mass parameter M^g slightly increases with 
the deformation parameters j3i of the two nuclei in the symmetric system ^^°Pd-|-^^°Pd as 
shown in Fig. 3. The small variation of M^^ is due to the shell structure. Therefore, the 
relatively large value of M^^ is a general result for collisions of both spherical and deformed 
nuclei. Fig. 4 shows the mass parameter M^^ at the touching configuration of symmetric 
systems = with /5j = 1 as a function of the mass number A. The mass Mg^ increases 
with the mass of the system. It slightly decreases with increasing mass asymmetry of the 
DNS. 



B. Diabatic two-center shell model 



Let us now consider the calculation of the mass parameters (pi]) with the diabatic single 
particle energies obtained with the method of maximum-symmetry in the diabatic two- 



center shell model p2| , |36| . In the diabatic motion the nucleons do not occupy the lowest 
single particle states as in the adiabatic case, but remain in their diabatic states. 

Due to the nonzero widths of the single particle states, the distribution of single particle 
strength over more complicated states is a Lorentzian distribution pfc(e) as introduced in 
Eq. (27) instead of the 5-function 5{e — e^) [|T8|. The occupation number n(efc) of a state k 
with energy and the corresponding value of fk are obtained, respectively, from functions 
n(e) and dn{e)/de calculated with zero widths of the levels 

n{ek) = / n{e)pk{e)de, (33) 



10 



A = -/^P.(e)cie. (34) 

The Lorentzian distribution increases the diffuseness of the Fermi-distribution. The Fermi 
distribution which is given at the touching configuration of the nuclei in the DNS is destroyed 
if the further motion of the system runs diabaticly. To treat the diabatic case, we use the 
following function n(e) for an arbitrary configuration of the system 

N 

nie) =Y.ai (e(e - q) - e(e - q+O) , (35) 

1=0 

where 0(2;) is the Heavyside's function and ei the energy of single particle state I with the 
occupation number a^. Here, the numbers / = 0, count the single particle states in the 
region of the Fermi level. The values Cq and ctv+i are the low and high limits of the single 
particle energies. For lower and higher energies, the occupation numbers are one and zero, 
respectively. Therefore, we assume ao = 1 and oat = in (^5l). The derivative dn{e)/de is 
expressed as 

dn(e) ^"^ 
^ = (1 - ai)S{e - ei) + ^ (a;_i - ai)S{e - ei) + aN-iS{e - cn)- (36) 

Then we obtain fk with (^^ as follows 

N~l 

/fc = (1 - ai)pfc(ei) + (a,_i - ai)pk{ei) + aAr-iPfc(eAr). (37) 

1=2 

In the calculations we assume the same average width for each Lorentzian pk{e). The diabatic 
occupation numbers ai are fixed at the touching configuration of the system using a Fermi 
distribution for a smaller temperature Tq* < Tq. The value of Tq is chosen in such a way that 
the occupation numbers n(efc) and the values of fk obtained at the touching configuration of 
the nuclei with the Fermi distribution or the approximate expressions (|35| ) - (|37|) are nearly 
equal. 

It is known that with zero width of single-particle levels the mass parameters are smaller 
in the diabatic case than in the adiabatic motion [18,22]. We found values for the masses 
Mxx and M^^ within the diabatic two-center shell model approach which are close to the ones 
obtained in adiabatic calculations because we take the width of single-particle states into 
consideration. The finite decay widths F^ destroy partially the diabatic motion (diabaticity). 
As result, the diabatic distribution of the single-particle occupation numbers approaches the 
adiabatic Fermi-distribution for the fixed values of e and A and excitation energy of system. 
The mass parameters at the touching configuration are the same in the adiabatic and diabatic 
descriptions because the single-particle levels and occupation numbers are practically the 
same in these two limits. As in the adiabatic two-center shell model, the mass parameters 
Mxx and M^e depend weakly on the mass asymmetry t] and on the deformations (3 of the 
nuclei [jS = jSi = (32 for ^^°Pd-|-^^''Pd) at the touching configuration. This result is presented 
in Fig. 5. Due to the regard of the width of the single-particle states, the mass parameters 
depend more smoothly on the mass asymmetry 7] than the mean field cranking masses. A 
stronger dependence of the mass parameters on rj is expected for entrance channels with 
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a larger mass asymmetry. In our previous work [16] we showed that the mass parameter 
of the neck could strongly increase with mass asymmetry starting from ri=0.7 even in the 
hydrodynamical treatment. In the present work we did not calculate the mass parameters 
for such very asymmetric systems because our version of the TCSM is not suitable for large 
1]. The mass parameter M^^ increases slightly with decreasing e or with increasing neck as 
shown in Fig. 6, is almost independent of the mass number A for heavy DNS and increases 
with A. In the reactions ^^"Pd+^^'^Pd and ^^Ca+^^^Hf, which produce the same compound 
nucleus ^^°U, the values of M^^ are similar for large e, but almost two times different for 
small values of e. 

The mass M^^ becomes smaller with growing temperature as shown in Fig. 7. The 
dependence of the mass parameter in neck on temperature Tq is caused mainly due to the 
width Ffc of the single-particle levels (F^ ~ Tq). The value of M^s would decrease as Tq^. 
The value of can be explained by the two-body contributions to the collective mass 

because Ffc accounts the nucleon-nucleon coUisions. So, Mf"^'(To=l MeV)/ M^'^f ^'(To=1.5 
MeV)~ 5 which is shown in Fig. 7. One- and two-body contributions [29] are contained in 
the nondiagonal and diagonal parts of the mass parameter M^e, respectively, are obtained in 
the same way. The one-body nondiagonal contribution to the mass is relatively insensitive to 
the temperature of system whereas the diagonal two-body contribution decreases strongly 
with temperature. However, for the excitation energy of the DNS formed in heavy ion 
collisions at low energy (< 15 MeV/nucleon) the main contribution to the mass parameter 
of the neck is due to two-body components. 

In Fig. 5 one can observe that the global minimum of M\x as a function of (3 in the 
reaction ^^°Pd-|-^^°Pd is around the ground state deformation of ^^"Pd (/? ~ 1.2). The mass 
Mxx grows with decreasing elongation A on average as shown in Fig. 6. The dependence of 
Mxx on A has fluctuations around an average upwards trend which are more pronounced with 
increasing total mass number A of the system (Fig. 6) and at smaller temperatures. These 
fluctuations are related to an increasing number of crossings between the single particle 
levels with decreasing A |3^. The diabaticity destroys the initial Fermi-distribution of the 



occupation numbers at the touching configuration and gives rise to fluctuating values of fk 
with decreasing elongation A and neck parameter e of the system. For larger temperature, the 
average width F increases and the function becomes smoother. Finally we note that the 
mass parameters remain always much larger than those calculated with the hydrodynamical 
Werner- Wheeler approximation. 



IV. SUMMARY 

Mass parameters for the relevant collective variables of strongly deformed systems con- 
sisting of two nuclei as the DNS are evaluated in the two-center shell model. Formulas for the 
masses are derived within the linear response theory by taking the fluctuation-dissipation 
theorem and the width of single particle states into account, and by other methods as Fermi's 
golden rule and a spectral smoothing in the mean-field cranking formula. The obtained mass 
parameter for the neck degree of freedom is much larger than the one obtained in the hydro- 
dynamical model with the Werner- Wheeler approximation. By applying the microscopical 
mass parameters we find a relatively stable neck during the time of reaction. In addition a 
large energy threshold due to the structural forbiddenness effect [0,0 hinders the motion 
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to proceed to smaller internuclear distances. We observe a strong correlation between a large 
mass inertia and a large degree of the structural forbiddenness for the heavy DNS as a func- 
tion of the fragmentation. Therefore, the DNS configuration exists a sufficiently long time 
and thermal fluctuations in the mass asymmetry coordinate could develop and lead to fusion 
and quasifission which are the essential reactions in the DNS concept. The comparison of 
our predicted results in the DNS concept with many experimental data for fusion reactions 
let us conclude that the DNS is stable against a melting in A and e. Our results also support 
the existence of long living cluster configurations which appear in ternary fission [p7|. 
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Fig.l. G.G.Adamian et al. 

FIG. 1. Dependence of the mass parameter M^^ (lower part) and of the shell correction SU 
(upper part) on e for the system ^^°Pd+^^°Pd at A = 1.6. In the calculation of M^^ an excitation 

energy of 30 MeV of the DNS and adiabatic single particle states are used. The mass parameter 
calculated in the Werner- Wheeler approximation is presented by a dashed line in the lower 
part. Units: m fm^ with m =nucleon mass. 



16 




10 15 

t(10"''s) 



20 




1.60 



1.62 



1.64 



1.66 



1.68 



1.70 



1.72 



Fig.2. G.G.Adamian et al. 

FIG. 2. Upper part: Time-dependence of the neck parameter e in the system ^^Zr+^^Zr cal- 
culated with microscopical (solid curve) and Werner-Wheeler (dashed curve) mass parameters. 
Lower part: Trajectories in the (A,£)-plane calculated for the system ^^^Xe-|-^^^Xe with micro- 
scopical mass parameters (solid curve) and with Werner-Wheeler mass parameters (dashed curve) . 
The end points of the solid and dashed curves in the drawing are at time t = 2 x 10~^^s and 
i = 2 X 10~^^s, respectively. 
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Fig. 3. G.G.Adamian et al. 

FIG. 3. Mass parameter M^g as a function of deformation /3 calculated for the system 
^^'^Pb+^^'^Pb at the touching configuration with excitation energy 30 MeV and adiabatic single 
particle states. Units: m fm^ with m =nucleon mass. 
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Fig.4. G.G.Adamian et al. 

FIG. 4. Mass parameter M^e as a function of the mass number A of symmetric systems {r] = 0) 
calculated with adiabatic single particle states for /?i = 1, A = 1.6, and £ = 0.75. Units: m fm^ 
with m =nucleon mass. 
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Fig.5. G.G.Adamian et al. 

FIG. 5. Upper part: Mass parameters M\x (solid line) and M^^ (dotted line) as functions 
of mass asymmetry r] at the touching configuration for DNS which leads to the same compound 
nucleus ^^'^U. The value of the neck coordinate is £=0.75 and excitation energy of the DNS is 30 
MeV. The calculation is done with spherical nuclei. Lower part: The same as in the upper part, 
but as a function of deformation j3 for the reaction ^^°Pd+^^°Pd (77 = 0). Diabatic single particle 
states are used. Units: m fm^ with m =nucleon mass. 
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Fig.6. G.G.Adamian et al. 
FIG. 6. Upper part: Mass parameter M^^ as a function of the neck coordinate £ at the touching 
configuration in the reactions: s^Zr+s^Zr (sohd hne), i^°Pd+^^°Pd (dashed hne), ^^exe+^^Sxe 
(dotted hne) and ^^Ca+^^^Hf (dashed-dotted hne). The excitation energy of the DNS in these 
reactions is 30 McV. Lower part: The same as in the upper part, but for Mxx as a function of the 
elongation A. Diabatic single particle states are used. Units: m fm^ with m =nucleon mass. 
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Fig.7. G.G.Adamian et al. 

FIG. 7. Mass parameter M^^ as a function of e at the touching configuration in the 
^^'^Pd+^^'^Pd reaction for temperatures Tq = 1 MeV (sohd Une) and 1.5 MeV (dashed hne). Dia- 
batic single particle states are used. Units: m fm^ with m =nucleon mass. 
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